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Semi-classical quantisation of space-times with apparent horizons. 
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Coherent or semiclassical states in canonical quantum gravity describe the classical 
Schwarzschild space-time. By tracing over the coherent state wavefunction inside the 
horizon, a density matrix is derived. Bekenstein-Hawking entropy is obtained from 
the density matrix, modulo the Immirzi parameter. The expectation value of the 
area and curvature operator is evaluated in these states. The behaviour near the 
singularity of the curvature operator shows that the singularity is resolved. We then 
generalise the results to space-times with spherically symmetric apparent horizons. 



I. INTRODUCTION 

Classical Black Holes, are observational realities, however the semi classical physics as- 
sociated with them remains to be explained. The black hole horizon is attributed with 
entropy, temperature, and a non-unitary form of thermal radiation or Hawking radiation [1] . 
What is the microscopic origin of entropy? How does a quantum mechanical wavefunction 
describe the horizon? Is the quantum black hole a pure state in the quantum theory? Why 
does the horizon radiate, and what is the end point of evaporation? Despite many plausible 
explanations, little has been achieved to demonstrate the complete truth. 

Further, classical general relativity predicts destruction of all material falling inside the 
horizon due to the presence of a central singularity. Do all the material accreting into 
the black hole perish at the central singularity even in a quantum mechanical description? 
The answers to these questions requires a full understanding of quantum geometry, as at 
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distances of the order of Planck length near the singularity, quantum effects will dominate 
over the classical prediction of a curvature singularity. 

To obtain such a quantum description, one needs a theory of quantum gravity. However, 
a complete theory of quantum gravity does not exist, though glimpses of truth have emerged 
in certain regimes. One such regime has been semiclassical gravity, where previously, quan- 
tum fields in curved space-time [1] were studied. Gravity remained classical. However, 
with the development of non-perturbative quantum gravity, relevant questions where one 
could 'semiclassically quantise' a given space-time have been answered to a certain extent. 
Semiclassical states have been constructed in canonical quantum gravity and we discuss the 
coherent states in this paper. These states have been well known in quantum mechanics and 
are 'wave packets' as opposed to exact eigenstates. They provide the closest approximation 
to classical physics as uncertainty is minimum here and the states are peaked in both the 
momentum and position representations. Thus expectation value of both momentum and 
configuration space variables are closest to their classical values as measured in these states. 
In case of canonical gravity, where time is separated, and the intrinsic metric of constant 
time slices and the extrinsic curvature constitute the phase space, these states can be used 
to build a entire space-time. Here a given black hole space-time with spatial slicings which 
include the horizon and the central singularity are discussed. To locate the apparent horizon 
and the central singularity one has to measure both the intrinsic metric and the extrinsic 
curvature of the spatial slice. This 'simultaneous' measurement is possible in the coherent 
states, as they are peaked in both the phase space variables, and the uncertainty is mini- 
mum. This would have been impossible in any other semi-classical state, where measuring 
the extrinsic curvature would have resulted in a complete uncertainty of the intrinsic metric. 

Coherent states for gravity were constructed in [7] using a formalism due to Hall. In 
this article, we address coherent states for black holes, first introduced in [11], and then 
used to derive the entropy of the black hole apparent horizon [12]. We give a complete 
derivation of a density matrix for the apparent horizon here by tracing over the wavefunction 
inside the horizon. The entropy is then obtained by the definition S = —Tip In p. This 
gives entropy to be proportional to area of the horizon in the first approximation, modulo 
a constant, which can be fixed to 1/4 due the the Immirzi parameter ambiguity in the 
formulation of the theory. The curvature operator or the Kretschmann scalar expectation 
value is also studied in detail in the states. The central singularity in the classical curvature 
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is clearly resolved in the semiclassical expectation value of the operator, mainly due to the 
uncertainty which prevents any measurement of area 0. Thus a upper bound exists for 
the curvature operator value proportional to the semiclassical parameter t which measures 
quantum fluctuations around a given classical geometry. As the semiclassical parameter 
goes to zero, the singularity re-appears indicating classical physics. 

The area operator is also examined here, as this crucially determines the entropy-area 
law. The classical value of the measured is given by equispaced numbers, and 

the entropy is actually proportional to the degeneracy of the area operator. The apparent 
horizon equation, which introduces correlations in the coherent state wavefunction does not 
impose any additional constraints on the area eigenvalues, and hence the counting yields a 
different value of the Immirzi parameter as obtained in [18]. 

In the first section, we review the coherent states, and in the next section we give an in- 
troduction to the classical phase space for gravity and discuss applications of coherent states 
to the same. The canonical variables, the black hole phase space, the corresponding coher- 
ent state, and a evaluation of the expectation value of the curvature operator are discussed 
next. The apparent horizon equation is examined in details, and a method of isolating the 
boundary conditions to be imposed on the coherent state wavefunction is analysed. The 
apparent horizon is a difference equation in the canonical discretised variables, and intro- 
duces correlations across the horizon. When the wavefunction inside the horizon is traced 
the density matrix describing the black hole space-time is described. The entropy and the 
Immirzi parameter are discussed, and the article concludes with a discussion and projects 
for future. The entire formulation here for the derivation of the entropy can be extended to 
include space-times with spherically symmetric apparent horizons. 

II. COHERENT STATES 

The coherent states are constructed to obtain classical physics from quantum mechanics. 
The origin of these states is well known in quantum mechanics for the simple harmonic 
oscillator, where the states appear as eigenstates of the annihilation operator a\z >= z\z >, 
(z is a label and represents a point in the complexified classical phase space x c \ — ip c \, x c \ 
denotes position and p c i are momentum). In general, according to Klauder [2, 3], coherent 
states are labelled by a continuous parameter z, and provide a overcomplete basis for the 
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Hilbert space, for a appropriate measure dfi(z). 

J \z >< z\dn(z) = 1 (1) 

A further restriction on the z to label points in classical phase space for a given system, 
fixes the state uniquely. For the Harmonic Oscillator, the Coherent states are also minimum 
uncertainty states or AxAp — h/2 (x,p being the configuration and momentum variable). 
Before going to the generalisation for gravity, I quote from [3] : Classical dynamics is quantum 
dynamics restricted to the only quantum degrees of freedom that may possibly be varied at a 
macroscopic level, namely, the mean position and the mean momentum (or velocity). This 
statement follows from the assumption that the classical action principle can be derived 
from a 'quantum action principle' if the states of Hilbert space are restricted to the coherent 
states. The quantum action principle is [2] 



' quantum 



< ip\i^-\ip > - < ip\H\ip > 
dr 



dr (2) 



\ip > is a wavefunction in the Hilbert space, r is the time parameter, and H the Hamiltonian 
of the system. The Schrodinger equation results with the variation of < ip\. However, if one 
takes the wavefunction to be of the form 

\,p >= e- iqp e ip Q\Q >, (3) 

where q, p are macroscopic position and momentum labels, and P, Q are corresponding 
quantum operators, then the above action principle, with the variation of p, q (hence the 
coherent states) lead to the following equation: 



-^resquantum 



/ [pq- <H>]dr (4) 



Clearly, if < H > is the classical Hamiltonian, then one recovers the classical Hamiltonian 
equations of motion. The observation of [3] is that the coherent states comprise a restricted 
set, and in all the exactly solvable systems, the states yield the classical system. In case of 
gravity, a different set of definitions have been used to generalise the coherent state structure. 
The basic features however remain the same. 

1) The states are labelled by points in the classical phase space, 

2) They provide a resolution of unity. 

The new definition of the coherent state generalised to gravity is due to Hall, [4]. In a 
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coherent state transform states in the configuration space (L 2 (R))aie taken to states in the 
Hilbert space defined on the holomorphic sector of the complexified phase space (H(C) R 
L 2 (C)). The kernel of the transformation is a coherent state wavefunction. The wavefunction 
is also the analytic continuation of the heat kernel of the Laplacian, which after appropriate 
normalisation corresponds to the kernel that appears in the kernel of the Coherent state 
transform [4, 5]. The coherent state obtained thus, coincides with the harmonic oscillator 
coherent state wavefunction, which is a eigenstate of the annihilation operator. The coherent 
state transform can be now defined for Hilbert spaces for arbitrary gauge compact gauge 
groups, in particular for SU(2). Canonical gravity in the Ashtekar-Barbero-Immirzi variables 
takes the form of a SU(2) gauge theory with additional constraints due to diffeomorphism 
invariance. The generalisation of the coherent states to a diffeomorphism invariance context 
appears in [6], and a complete study of their properties are in a series of papers in [7, 8, 10]. 
A very interesting question arises: Is there a similar 'Quantum Action' principle as in (2) 
for Gravity, and if so do the coherent states restrict the action to the 'Classical action' as in 
(4)? This indeed is a very difficult question, as no one knows a 'corresponding' Schrodinger 
equation for the quantum gravity states, and the only known 'equivalent' equation, the 
Wheeler- Dewitt, equation is derived from the classical action, (5S/5N = 0), where N is the 
lapse. Hence would not qualify as a ab initio 'quantum equation' (The Schrodinger equation 
is not derived from the classical action for quantum mechanics). For the Hall coherent states, 
the SHO coherent states are recovered in case of L 2 (R), and Equation (3) is indeed true. 
However, for the coherent states of gravity, which are defined on SU(2) group, what would 
be the corresponding equations for (2) and (4). We discuss this in the next section, after 
identifying the phase space for gravity. For the sake of clarity, we define the coherent state 
as [4, 5], 

^\z)=p\x) x ^ z (5) 

where, p is the heat kernel for the Laplacian on the given configuration space, (for SHO the 
configuration space is R). z takes values in the complexified phase space, and corresponds 
to the continuous label of states as per the definition in [2], and the other parameter t, 
is essentially the 'semiclassicallity' parameter. This parameter gives the width or variance 
of the coherent state around the mean value or peak value in position space. In case of 
the simple harmonic oscillator, expectation values of operators correspond to exact classical 
values, and hence, < P,Q >= p,q, irrespective of the semiclassicality parameter. For the 
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coherent states in gravity, these statements are true only in the limit t — > 0. Hence, any 
expansion of the expectation value of operators in this parameter t is actually a study of 
quantum fluctuations around a given classical geometry. The parameter t, as defined in [7] 
is defined as ti/a, where l p is the Planck length, and a a dimensionful parameter, which in 
the case of Schwarzschild black hole can be r 2 g (horizon radius squared). 

Continuing the discussion on the coherent states, the generalisation of the above definition 
in Equation (5), leads to the following for the Hilbert space of any arbitrary gauge group 
H, whose complexified elements lie in G. 

^ t {g) = p t {h) h ^ g (6) 

= exp(—£V)5hh' is the heat kernel of the Laplacian on the group manifold.) with an 
appropriate normalisation (h C H,g C G). The states are overcomplete with respect to a 
measure dfi(g) which in the case of SU(2), was shown to be the Liouville measure [4, 7]. The 
Laplacian for SU(2) corresponds to the Casimir operator which has the eigenvalues j(j + 1) 
in the jth irreducible representation, the coherent state can be written as a sum over the 
irreducible representations, using a theorem due to Peter and Weyl. 

i>\g) = Y. d ^ m+1)/2 xA9h- X ) (?) 

3 

(dj is the degeneracy of the irreducible representation with character Xj)- Since it is this 
form which is relevant for Canonical gravity, we proceed to find appropriate phase space 
variables for gravity, and then define the coherent state as a function of the phase space 
variables. 

III. CLASSICAL PHASE SPACE FOR GRAVITY 

We study gravity with the space-time metric g^ as the configuration space variable. 
Due to diffeomorphism invariance, the reduced space is really g^j Dif f(M). A separation 
of time and space in the ADM formulation further fixes 'the configuration space' as the 
intrinsic metric q a b, of the constant time slices, the lapse and the shift for propagation along 
the time like directions, are given by N, N a , where a = 1, 2, 3. One can define the momenta 
conjugate to these variables from the classical action. As is well known, the Hamiltonian, 
which is dual to the lapse N, is a constraint in gravity due to the absence of the N term in the 
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action. Thus, only on the constrained surface, the actual variables are q ab and the canonical 
conjugate variable 7r afe = q' x l 2 [K ah — q ah K') , where K ab is the extrinsic curvature of the slice. 
A coherent state for these geometrodynamical variables is yet to be constructed. To use the 
Hall coherent state, one has to use the new variables formulation of canonical gravity. Here 
the tangent space of each point on the spatial slice is used to define the variables: 

A{ = T T a - (3K ab E bI , PE^E'j = detq q ab (8) 

with (3 being the arbitrary parameter in the theory or the Immirzi parameter, and I runs 
from 1,2,3 to denote the SO (3) or SU(2) degrees of freedom of the tangent space. is the 
spin connection, E bI are densitised triads and A ! a is the SU(2) gauge connection. 
In these, the action for gravity has the form 

I G = - [ d z x{^} - A 7 G/ -NH- N a H a ]dr (9) 

The above has the form of a Yang-Mill's action, however with additional constraints in 
the form of the Hamiltonian H, and the diffeomorphism generators H a (G 1 is the usual 
SU(2) Gauss's constraint, A 7 is the Lagrange multiplier). On the classical phase space, 
the constraint equations are H, H a = and the canonically conjugate variables are A ! a , Ef. 
Clearly no one knows what the action principle for Quantum Gravity is (or the Schrodinger 
Equation) is for quantum gravity. Does the restriction to Coherent states yield the classical 
equations for the canonical variables? As we know, that finally the quantisation of the action 
in Equation (9) is carried out in the smeared variables called the holonomy: h e (A) which 
are path ordered exponentials of the gauge connection along one dimensional analytic edges 
e, and the corresponding dual momentum P\ which are the densitised triads smeared along 
2-dimensional surfaces. This is effectively a discretisation of space, and one defines the basic 
variables over graphs, and their duals. In [7], Thiemann explored the variables originally 
defined by Ashtekar and Lewandowski further, and succeeded in 'quantising' these graph 
dependent variables, and also obtained the appropriate classical limit for these operators. 
The variables are: 

h e (A) =Vexp([ A) Pi(A,E) = -Tr[T T h e f h^Eh^h; 1 } (10) 

J e (X J S 

(T 1 = — 27 1 are the generators of SU(2)(7 / are Pauli Matrices), and a is a dimensionfull 
parameter, usually fixed as a function of the parameters in the classical theory). The h p 
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are the holonomies along edges defined on the 2-surface. The variables satisfy the Poisson 
Algebra, 

{h e , h e ,} = 0, {h e , P$} = -Jee'h^- {p/„ P/} = ~/ 3K P? $ee> (11) 

The complexified element formed from these phase space variables (similar to x — ip in R) 
is an element of SL(2,C), and is g e = e~ lT ' p ' ^ 2 h e . The Coherent state is constructed to be 
peaked at the SL(2,C) valued element. The classical action for these discrete variables (for 
a particular edge) will be (on the Constrained Surface): 

s = Z J XI (12) 

If we define the 'quantum action' principle for gravity in a similar way, except that one 
confines oneself to the derivation of the Kinetic term (12), as above, then, it has the form: 

S q = ij 'dr<i\^\i> (13) 

where < £| is a arbitrary state in the Quantum Hilbert Space, then does one recover (12) 
by confining oneself to the coherent state? Or in other words, what is 

f d 

Sresquantum = % J dr < if)* | — > (14) 

(where i/j* is the coherent state for a single edge.) Strangely enough, though in a rather 
straightforward calculation (reported in the Appendix(A))) , one recovers (12), in the limit 
t -> 0. 

So in principle, we have succeeded in deriving a appropriate phase space for gravity, 
for which a coherent state can be defined, in a very similar manner to any other quantum 
mechanical system, and they give classical physics. All the previous discussion is about 
the coherent state peaked at classical phase space variables of one particular edge. The 
entire manifold is however charted with a graph, comprised of edges linked at vertices. A 
SU(2) Hilbert space is associated with one edge, and the complete description of the entire 
manifold, is a tensor product of Hilbert spaces of all the edges comprising the graph. The 
coherent state for the graph will thus be of the form 



e 



(15) 



9 



This is a gauge covariant tensor product, of the Coherent state defined on each edge. (The 
gauge transformations act on the holonomy and the corresponding momenta thus: 

h e - g(0)h e g{l)-\ Pi - g^PigiOy 1 (16) 

The g(0) and g(l) are the SU(2) valued group elements acting at the starting point and 
end point of a edge respectively.) The gauge invariant coherent state has intertwiners at the 
vertices, which ensure that the state transforms as a singlet at the vertices. However, for 
this article, we will confine the discussion to the gauge covariant coherent state, and in the 
conclusion, comment on the complications which can arise in a gauge invariant state. 

To end the discussion of the coherent state for gravity, we emphasise the following two 
labels: 

1) The classical phase space label g e . These variables satisfy all the constraints by construc- 
tion. Despite the discretisation involved through the definition of a graph, the variables 
respect the inherent continuity in the classical metric. 

2) The semi-classicality parameter t. 

The actual Coherent state is defined over the tensor product of Hilbert spaces, which for 
asymptotically flat manifolds can require a infinite tensor product of Hilbert spaces [10]. 

IV. SEMICLASSICAL BLACK HOLES 

A. The classical phase space 

The classical phase space for the spherically symmetric sector in gravity, will be consti- 
tuted by Schwarzschild black holes. Clearly, one already knows the spherically symmetric 
classical solution of Einstein's equation in vacuum, and for this metric, all the constraints 
are satisfied. The intriguing part is to isolate the graph dependence to derive the discrete 
classical phase space of Equation (10). What exactly would be an appropriate graph? Ob- 
viously whatever the graph, the classical discretisation has to be spherically symmetric. 
One very convenient set of graphs is to take edges along the coordinate lines r, 9, <fi with 
appropriate discrete labels attached to them. This was what was done in [11]. The dual 
polyhedronal decomposition is then comprised of spherical surfaces which intersect the edges 
at their middle points. The holonomy and the momenta are calculated in [11], and their 
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behaviour analysed. The radial holonomy and the corresponding radial momentum have the 
expressions: 



h f 



1 



cos r 



1/2 

T 2 



,1/2 



ry 1 sin t' 



1 



1/2 

r 2 



,1/2 



(17) 



Where T\ is the begining of the edge, and r 2 is the end of the edge, (r' = y^/2 for /3 = 1). 
The holonomy is ofcourse independent of the angular coordinates of the begining and the 
end of the edge. Though the holonomy depends on the extrinsic curvature which diverges 
at the singularity, there is no such divergence in the regulated 'holonomy'. 

'e(l) 




The evaluation of the momenta is a little complicated, as it involves the evaluation of 
the integral on a two surface, with the triad convoluted with the holonomies associated with 
paths from a generic point to the point at which the edge intersects the two surface. As 
given in the above figure, the point of intersection of the edge with the surface is denoted as 
O and has the coordinates (r,#o,0o) ; {i n the final expressions 0o does not contribute}. The 
width of the surface is from 6 — 6,6 + 8, and a linear dependence on the width in the 
direction, which is suppressed here for brevity. The width along the radial edge is given by 
5. For details see [11], 



P 1 



-Tr 



7 ' ( cos vi) ~ %lX sin vi) ) x(r) ( cos vi) + ll1 sin 



2r 



X(r)=X 1 (r) 1 1 + X 3 (r 



= 7 



X, 



a 



with 



X 1 (r) 
Mr) 



r 

—7- sin 9 { 

or 



r 

—7 cos 6 

or 



sin(l - a')6' sin(l + a')9'~ 
l-a> + (1 + a') 

sin(l - ay _ sin(l + a')9 r 
(1 - a') (1 + a') 
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(18) 
(19) 

(20) 
(21) 



Where a 1 



a 2 +l 



2 ' 



a 



—. Now, one uses the above in (18), to get the following result 



for the momentum components 
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Xr r / i a\ ■ ( i • 

-a cos 1 7 a J + sin (7 a' 



( 7 ' = 5/2r 5 ). 



^V^ + ^ + ^f = I^XJ+XI (23) 
9 



a a 4 



/ sin [(l-a')0]\ 2 /sin[(l + a')0]' 
I (1-a') J + { (l + «0 , 



- 2cos(2. ) Sm ^-^ S1 ^ 1 + ^ nl/2 
1 — a' 1 + a' 

Also, this complicated dependence on the coordinate point 9q at which the edge intersects 
the dual surface becomes clear, when one takes the size 9 of the surface to be very fine. One 

notes that as the graphs get finer, the above approximates 9 — > (restoring the width 0): 

2 

P er = 2-sin0 o 00. (24) 

(X 

Which is the area of the 2-surface which a edge intersects at it's middlepoint. Clearly, 
the classical value of area is given as above. Question is, when one lifts the variable 



P/P/ to an operator, is it diagonalised in a area eigenstate? The answer is yes, in a 
exact orthonormal eigenstate, this corresponds to the area operator and has the eigenvalue 



j(j + \)t. However, when we evaluate the expectation value in the coherent state, classical 
area is corrected from this exact eigenvalue, and in the limit t — > this has a different 
spectrum as we observe in Section (IV C). Before one discusses the coherent state defined 
for these variables, some crucial points need to be noted: 

1) The Holonomy and the Momenta remain finite numbers even in the vicinity of the 
singularity. This can be attributed to a regularisation achieved due to the discretisation. 

2) The variables are continuous across the horizon, and contain information about the 
apparent horizon. 

3) For very fine graphs, one obtains as P the classical area of the bit of the dual surface, 
induced by the classical metric. 



B. The Coherent State 



The coherent state for the black hole phase space can be written explicitly as in [8]. 
The states are peaked at the classical values, and are well behaved in the entire black hole 
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slice. Due to the fact that the classical variables h e , P e are themselves well behaved, there is 
no large fluctuations in the quantum corrections, at the horizon, or even at the vicinity of 
the singularity. This is a indication of the 'minimum uncertainty' principle obeyed by the 
states. The deviation around the mean value of both the momentum, and the holonomy are 
measured by the variance, and they are such that 

T i 

< t/j^Ahe AP e V >= — h e (25) 

which being the classical holonomy, by Equation (17) is always within the value —1... + 1 
in magnitude. Hence, the Coherent states, are ideal to study the semi-classical regime, and 
due to the 'regularisation' achieved in the set of variables, here, the semiclassical approxi- 
mation is valid in the entire black-hole slice. The coherent state in the configuration space 
representation has the following expression 

^{ghr 1 ) = E(2j + ^e-^+^igh- 1 ) (26) 

3 

where j labels the eigenvalue of the SU(2) casimir, and Xj corresponds to the character of 
the corresponding irreducible representation. The corresponding momentum representation, 
is determined by taking a 'Fourier transform' in the character space of the irreducible states, 
and one obtaines 

< tf\jmn >= e-^+^TT^W (27) 

These are the position and momentum representations of the coherent state peaked at in- 
dividual edges, with the classical values encoded in the SL(2,C) valued variable g e . In this 
paper we make a very interesting observation for the state which is 'peaked' at the classical 
value of P = or area =0, with arbitrary holonomy. For this state, the expectation value of 
the area operator is however non-zero, and proportional to t, as expected from the minimum 
uncertainty criteria demonstrated in equation (25). 

C. The Equispaced Area spectrum 

There appears to be two different questions regarding the spectrum of the area operator. 

1) The eigenvalue of the area operator as obtained in an exact eigenstate. 

2) The spectrum of the area of the black hole horizon, as measured in a appropriate quantum 
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state or a semiclassical state. 



Regarding the first question, the usual regularisation of the area operator [15] gives the 



eigenvalue of the area operator intersected by a edge to be : 8ir^Jj(j + 1)1^ in a kinematical 
eigenstate of the same operator [15]. In this context, in [13] it is claimed that the SU(2) 
casimir, whose square root is proportional to the area operator undergoes a renormalisation 
and gives the area eigenvalue as (j + 1/2)1^, in it's eigenstate [13]. a Previous computations 
of entropy , counted the degrees of freedom of a boundary Chern Simons theory at the 
horizon, given that the area of the horizon assumed certain values as measured in exact 
eigenstates associated with edges crossing the horizon. These edges as consistent with the 
previous calculations induced the horizon with bits of area SftyJij + 1)^. 

The calculation which explains the semi-classical processes like black hole entropy or 
Hawking radiation should arise from a appropriate semi-classical limit of a quantum theory. 
To recover 'semi-classical' entropy one needs to find the microscopic degrees of freedom 
corresponding to a given classical black hole space-time. The only well known states in a 
quantum theory which give classical physics are the coherent states. So, in a coherent state, 
the classical horizon area should be the 'expectation value' of an area operator. As discussed 
in Section(IV A), the variable P corresponds to classical area. It is interesting to study the 
momentum representation of the coherent state wavefunction, which is expanded in the area 
eigenstates to determine the classical clXGcl clS cL function of the area eigenvalue or the SU(2) 
Casimir eigenvalue j. 

\i> >=T, d ^~ tKj+1)/2 ^(9)m.n\jmn >, (28) 

3 

where iij is the jth irreducible representation of g which is an SU(2) valued matrix encoding 
information about the classical variables. The state \jmn > is an exact area operator 
eigenstate. The coherent state is a superposition of such eigenstates with the coefficient 
being the spin j irreducible representation of g. Now, the probability distribution is, 

g 3 3+ ) |7Tj(^) mn | , 1 | 1 / ^ + D ^2 



ocexp(--((j + l/2)t-P) 2 ), (29) 
where P is the classical area. Hence, the classical area is P = (j c i + l/2)t, with maximum 



probability. It is also the expectation value of the operator y A + i 2 /4 in an exact eigenstate, 
where A is the usual area operator. Thus the classical area has a 'corrected spectrum' 
determined by the equispaced discrete numbers (j c i + l/2)t [11, 12]. Note that the word 
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'spectrum' has to be used with caution as this is actually the expectation value of the 
area operator in a coherent state. More precisely, given that the area operator is A, the 
expectation value of the operator in the coherent state is of the form: 

< A >=< ^\A\^ >= pp E yjj(j + i)t (2j + D^di) < !rj - u ~An) J (30) 

Using g e = h e e- iTlpI / 2 and |^|| 2 = ^f^V^/ 4 and ^(gfg) = Xj (H*), H = 
e -iT I P'/2^ _ gjjjjj^j -f 1)P/ sinhP. the following can be derived 



<A> = 3i nh P J (2^ _ * (2 . + 1)e .„, )(2J+2)/4 E ni{g1)M 



mn 



•£ ( 2J + 1)* exp (JBi + PV-^ 5inh < 2 t + " P 

2Pv^F ^V4 v ' y ; 4 v J ; P V 4t / sinhP 

g v^3^ n exp ( J n - 2P ^ \ ; n =(2j + l)t 



1 /" ? / (,x-2P) 2 \ , 

— ==- \ v £ —tx exp ax 



4P(2 

= 4P~ / ^ x ^ X ~ 2P ) rfx ( Limit * °) 

= P (31) 

Now, when the sum is converted to a integral in a variable x (step 3 of above), the corrections 
are proportional to t, and hence in first order in t, this is a perfectly valid result [21] 
. However, for even a detectable finite t, this result is true iff 2P = (2j c i + l/2)t, or 
P = {id + l/2)i, where j c i are discrete numbers. This observation is not contradictory 
to the previous work on area spectrum of loop quantum gravity, where the measurement 
is in an exact eigenstate. Note, in the t — > limit, to assign a finite area, one has to 
take jd » 1, and hence the equispaced spectrum assigns a continuous value as the non- 
equispaced spectrum : j d t to the area. 

When one is trying to count the number of ways to build a macroscopic area, using 
the coherent states, the classical areas induced by the edges, are counted by the discrete 
numbers jd- There is a degeneracy associated with every jd, corresponding to the numbers 
m,n— —jd--jci (corresponding to the expectation values of the operators P|,P|), and this 
is 2j c i + 1. This degeneracy gives a degeneracy associated with a given horizon area. Why 
this should be the entropy of a black hole is the subject of discussion in the next few sections. 
In the next subsection however, we discuss the resolution of the singularity at the center 
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of the black hole. We will now concentrate on the specific coherent state, whose classical 
label has contributions only from the unitary component of g e , namely the holonomy, and 
the P 1 = 0, or the classical area induced by the particular edge is 0. (or g e = h c i). In this 
particular choice of graph, this will happen for areas induced by edges close to r = 0. The 
coherent state is: 

l^oM >= E dje-W+WirjihaUnljmn > (32) 

jmn 

In the above g e = h d , as the P 1 = 0. The expectation value of the area operator in the 
coherent state is 

<^{h d )\P\^{h d )> 

= < ^o(^)l I! ^VjCTTTyte-^^+^^^^ljm^ > (33) 

jmn 



jmn 
1 



E*0' + l)%" JlJ+1)1 %(l) (35) 



2 

3 

= \t + 0(t") (36) 

In equation (34), one uses < jmn\klq >= j- S jk 5 m i5 nq , and TTj(h d ) mn = Hj{h ] d ) nm 

Any measurement on the coherent state gives the classical expectation value only when 
t — > 0, even a tiny amount of t would ensure that there is minimum area which one can 
measure in the coherent states. Thus one obtaines in some sense a minimum radius which 
can be measured in the spherically symmetric coordinates, 

< P>=-r 2 min sin e<fl(j> (37) 

CL 

r 2 

for each edge. Which gives the value of -°^ L = jp-t, where n is the number of edges inducing 
the total area of the sphere. Thus even for a infinitesimal t, and fine graph, there is a 
minimum radius and this uncertainty leads to a resolution of the singularity as we shall 
subsequently observe. 
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D. The Curvature operator and the Singularity 

The Curvature scalar R^ v \ a will consist of contributions from the intrinsic metric, as well 
as contributions from the extrinsic curvature. These are of the form: 



Clearly, when the intrinsic metric is taken to be flat, one has the curvature scalar as 



(38) 



K* - K bc K ac K d a K h d 



(39) 



(In principle there can be order t corrections to the intrinsic curvature, but to first order, the 
contribution is essentially zero. Hence we ignore the intrinsic curvature terms here.) This 
quantity diverges at r = classically for the Schwarzschild metric. One can write the entire 
above expression in terms of the regularised holonomy and momenta. To see the form, we 
write the expression for the extrinsic curvature in the form: 

K ab = ^ [Ti - A a ] EwV&tiE 

Further, the spin connection is written in terms of the momenta: 

1 



(40) 



JJK pb 
j_ 1 JJK rpb 

+ ^ E K 



E tb 



E, 



b,a 



2E 



J (detE), b 



E j (detE), 



(41) 



detE detE 

The entire curvature is then finally written in terms of the two operators h e (A), P e {A) and 
their expectation values in the coherent states. The R 2 operator is 

2 



R 2 



J7>a rpb ~\ ' 



(ri - K) (r /6 - A Ib ) (r k c - A k c ) E d k (v l d - a\) e\ (e?e?) (det e 2 



(42) 



Before one actually lifts the above expression to a operator equation h e ,P e , one must also 
express inverse powers of triads which appear in the Equation(42), in terms of Poisson 
Brackets. Thus a measurement of R 2 in the coherent state will be obtained after one has 
replaced the A a and the E T a in terms of the holonomy and the corresponding momentum, 
and the regularised expression for the inverse triads Ef. In the vicinity of the singularity, 
the terms containing Y ! a donot contribute to the singularity of the curvature but all go to zero 
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at r = 0, and hence can be ignored in the calculation for the upper bound of the curvature 
operator. It is the terms containing which are potentially divergent. Thus in the curvature 
operator, we retain the terms independent of the spin connection. Note, as discussed later for 
the apparent horizon equation, it is possible to write the extrinsic curvature operator solely 
in terms of the gauge connection operator, and the triad operator, by using the Immirzi 
parameter (100). Hence all the derivations for the extrinsic curvature, here will be true, 
when one takes into considerations the appropriate (3 of the theory. However, for the next 
few discussions, we 'ignore' the spin connection in the computation of the curvature operator 
in the vicinity of the singularity. Since the quantum fluctuations are always small, one can 
never induce large values for the spin connection operator, when their classical value is 0. 
Thus for a measurement of the curvature operator in the coherent state, in the vicinity of 
the singularity, it is always justified to ignore the spin connection operator. 

Before evaluating the expectation value of the curvature operator, we make the following 
observations on the operator ordering ambiguity which occurs for the operators which are 
functions of both h e and P/. These will clarify some of our assumptions and calculations. 

In field theories, operator ordering ambiguities often lead to infinities. However, for these 
coherent states, the ambiguities are proportional to t and should go to zero in the classical 
limit. This observation might not be true when the classical geometry itself has a singularity. 
To investigate the situation where the classical geometry itself is singular, we take a arbitrary 
function of / = hP with, the normal ordering defined as 

<V| :/: \^> = <^|(P/)/^> + <#(P/)A]|^> 

= (PX)d + \ T\h e ) d + 0(t 2 ) (43) 

Clearly, the operator ordering ambiguity is proportional to h ec \. Let us examine the radial 
case, Clearly, since by Equation (17), we find that taking one of the end points of the edges 
to ri(r 2 ) — > 0, 



T 



T 



T 



r 



cos 



COS 



T 



i 



COS 



sm 



r 2 



sm 



— sm 



n 

t' 



sm 



cos 



r 2 , 
t' 



(44) 



the limits cos(l/wr 1 (r 2 )), sin(l/-i/r 1 (r 2 )) oscillate within the finite limits — 1..1. Hence in 
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terms of the 'regularised variables' taking the limit to the singularity does not affect the co- 
herent state, and as such, indeed the coherent state can be defined in terms of the regularised 
classical variables, even in the vicinity of the singularity. Thus now any arbitrary polynomial 
in terms of the holonomy and momenta, can be obtained as a expectation value and classical 
limit determined even at the singularity. Thus we proceed to write the curvature operator 
in terms of the holonomy and momentum, and find the expectation value. 

It is a interesting calculation to realise how the classical singularity can be recovered from 
a 'finite regulated' holonomy, For example, with Equation (17), the radial gauge field e.g. 
will be: 

A ' = -w^F r{K - 1)1 (45) 

Where r\ and r 2 are the begining and end points of the edge respectively. Now the most 
general form of the holonomy can be taken as h er = e Tl(jI = cos a + "^sL s j n <j ; which implies 



1 „ r „r a 1 e IJK a J T K 



A 1 = ; -Tr [cos aT 1 sincr , 

2(ri -r 2 ) cr a 

1 a 1 
= sin a 



sin a\ 



ri — r 2 cr 



(46) 

Here the interesting observation is that now in the particular value of holonomy calculated 
above: a 1 cx (r^ 2 — r 2 2 )5 n , which by itself is divergent in the limit r 1; r 2 — > 0. a 1 jo = 1, 
and the sincr factor infinitely oscillates from -1..1 as a is diverges. Thus, even in the value of 
A r , the divergence does not appear to show up, in the regularised variables or the extrinsic 
curvature of the manifold if there exists a minimum edge length. However, it is the question 
of taking limits, and if somehow, one shrinks the edgelength faster than the limit r — > 0, 
then the sin a factor can be approximated as a, and the following divergence appears: 

T g -1/2 



ri - r 2 



r 

(47) 



In the limit t\ — > r 2 and r 2 — > 0, this is indeed a singular value of the extrinsic/gauge con- 
nection recovered. Thus it is very clear that in terms of the regularised variables, even the 
existence of a minimum edge length gives a upper bound on the curvature proportional to 
jz, where 5 is the minimum edge length. As described above for a non-zero value of t, there 
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is a minimum area one can measure in the coherent states, which would also give a bound 
on the length of the edges as measured. However, we will attempt to measure the curvature 
directly using instead the regularised extrinsic curvature operator, and demonstrate a pos- 
sible resolution of the singularity in terms of the coherent states. First let us approximate 
K a b, (ignoring T^) using the usual regularisation as follows: 



Using the previous regularisations of the inverse triad, one can obtain a suitable expression 
as follows: Writing el = CTr[T T hjMh eb ^ V}}, note that C is a graph dependent constant, 
which is necessary, as the volume operator, defined in terms of the PjS is graph dependent. 
In fact as in [14], the constants are fixed here by calculating the Poisson bracket of the 



holonomy operator with the volume operator. Since V = J ■^e IJK e abc PfP 1 jP^- , the Poisson 

bracket of the holonomy with the Volume operator is: 



K ab 




(48) 






- zy{ he ' y eIJKe abcPlP b jPK} 

— £ ^^Enhrhp P^P^ 

8aV J K 

= K ^h e j\ptr 



here, v = 9(p5/a 3 ^ 2 and s er = 9(p/a,s eg = 6(f)/a,s e , = 59 /a. (5,9,4> denote the edge lengths 
along the coordinate directions). Thus, now multiplying by T T h~^ and taking Trace, gives 



the constant to be 



C, 



2as 



(49) 



VK 



In the quantum commutator, C oc l/ll. 



K ab = C^TrlT 1 



2(e(l)-e(0)) 



}Tr[T 



(50) 



<^\K ab \^ >=C eb <i,\Tr[T I 



2(e(l)-e(0)) 



]Tr[T I (V-h- l Vh eb )]]\^> 



(51) 
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However, to avoid the presence of the double trace in the operator and one instead uses the 
following regularisation of the extrinsic curvature: 



K, 



ab 



4V 



a 



-Tr 



he a K b l {he b ,V} 



(52) 
(53) 



2(e o (0) - e a (l)) 

(e a (0),e a (l) denote the begining and end point of an edge). With the classical value re- 
covered in the limit the edge length goes to zero. The h ea denotes the holonomy along the 
edge e a , and V is the corresponding volume operator. The Poisson brackets give the inverse 
triads and the limit the edge lengths go to zero, one is left with the classical expression for 
the extrinsic curvature. However, when one lifts the above to an operator equation, the 
operator ordering is taken to be: 



K, 



a 



ab 



'-b 



-Tr 



K'[h eb ,v]h ea 



(54) 



2(e o (0) - e a (l)) 

As this ensures that the diagonal components are recovered appropriately. For the radial 
component of the extrinsic curvature, one obtains: 



K r 



-Tr 



h- e }[h er) V}h er 



(55) 



2(e r (l)-e r (0)) 

Thus a typical term in the evaluation of the curvature will include the Volume operator. 
To find the spectrum of the volume operator in these coherent states, we have to derive the 
operator in more details. 



E. Volume Operator 

As in contrast to previous derivations of the volume operator, here, there exists a classical 
metric, to fix the constants and the exact expression for the operator in terms of the gauge 
invariant vector fields. We will proceed in the following manner: Firstly, we fix the classical 
volume in terms of the graph degrees of freedom, and then lift the obtained expression to 
a operator in the quantum theory. The volume operator will be obviously adapted to the 
specific graph chosen. Here we will also focus on a set of vertices in a local region R. The dual 
polyhedronal decomposition of the manifold will be important in the determination of the 
individual volume cells. Here, the graph has been taken to be 6-valent with three ingoing 
and three outgoing edges at a given vertex. Since the classical intrinsic metric has been 
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taken in the spherical coordinates, the edges are aligned along the coordinate directions. As 
described before the triads are smeared over the dual 2-surfaces which the edges intersect at 
their midpoints. The dual surfaces constitute a volume cell, with a vertex at the center of 
the cell. This geometric construction very conveniently follows from the definition [7]. The 
figure shows clearly the construction of the volume cell, and each vertex is surrounded by 
the dual surfaces forming the walls. The volume associated with each vertex is therefore 



(56) 



Where e„, triplet of edges intersecting at the vertex. For a generic vertex located 

at the point (as measured by the classical metric ) (r, 8$, 0o)> the classical volume as evaluated 
from above can be evaluated. Typically, For the edges ingoing at the vertex, contribute with 
(23) [ for the momentum of angular edges [11]], 



Pi = (r- 5f sin(£ - 0/2)^ P* g = (r - 5/2) sm(9 - 



(r-5/2) 6 - (57) 



The contribution from the outgoing edges are similarly of the form: 



Pi = (r + 5f sm(6 + 0/2)^ P* e = (r + 5/2) sm(6 + 6)^ P^ = (r + 5/2)^ (58) 



The determinant to first order in the edges of the cube (5, 9, <f>) is obtained as: 



V = r 2 sin^o 



7,3/2 



(59) 



Which is the required volume of the cell. 

e 

e 




e 

e 
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Thus the entire volume of the manifold can be obtained as the sum of the volume of 
the cells. The spectrum of the volume is then obtained in the coherent state by using the 
standard techniques. The operator P/ is replaced by the Right invariant vector fields Xf 
and then the expectation value of the operator is obtained in the state. This operator is of 
course quite similar to previous derivations, however, here the basic cell is prompted by the 
classical metric, and the six-valent graph adapted to that. Now, to evaluate explicit matrix 
elements, in the coherent states, one obtaines the following: 

< V 2 >= 1 < ^euKXlXlXfcW > . (60) 

This operator, again due to a trick due to Thiemann, can be written as: 

e IJK X?X*X c k = \{X a + X c ) 2 , (X b + X c ) 2 } (61) 

which however is not a convenient set of operators at this juncture, as we are interested 
mainly in the classical limit, where a naive evaluation of the operator expectation values 
gives 0, (recall that due to the nature of the classical metric, < X a .X b > =0). Instead, we 
use the sperical symmetry of the classical metric to equate the following: 

Limit^o < V >= Limit t _ < P r > — 7= (62) 

V a 

as P r = r 2 sin 9 9(f) has the same magnitude as Jq. This would be true only at the first order 
in the metric, and the result will be considerably different in higher order corrections in the 
semiclassicality parameter t. However, since we are interested in a plausible resolution of 
the singularity, we approximate the volume operator at this level of the discussion by the 
above. 



F. Resolution of Singularity 

To find the expectation value of the curvature operator in the coherent state, we evaluate 
the explicit expressions of the complete operator. The operator is taken as a density 1 
operator as this gives some nice properties. Thus, the operator is of the form: 

y/gR 2 (63) 
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by construction. This as by equation (39) now has the form (including only the potentially 
diverging terms): 



2^ 



K 4 - K ab K cd K ac K bd 



2^ (KaoK^qOy - K ab K cd K ef K gh q ae qf c q b °q° h 



(64) 
(65) 



The regularisation of the extrinsic curvature is then used as from equation (53), however, 
the inverse metric q ac , when regularised, introduces inverse powers of detg, which are then 
absorbed into the Poisson brackets in the numerator, and also with a point splitting method 
introduced by Thiemann. This essentially has the following method, which uses the fact 
that: 

detei 



e IJK e abc e a je b je c K 



1 

3!^ 

C a ChC c tjk abcrp 

3!^ 



Th^{h ea ,V} Tr T J h-'{h eb ,V) Tr T K h^{h ec: V] (66) 



The above Poisson brackets then add to the brackets with the volume, and the inverse 
powers in the denominator can then be absorbed in the numerators. We concentrate in the 
first term of Equation (65), and an explicit expression with the extrinsic curvatures is the 
following: 



. detel 



,1/4 



(67) 



Now, we find using the point splitting of equation (66), and the regularisation of the Extrinsic 
curvature of Equation(54): 



dete* 



a K 2 



q 

(~* (~* 



3!gV4 
1 



T J h^{h em , V}\ Tr [T J h^{h en ,V}\ Tr [T K h^, {h £p , V} 



x 



-Tr[h-Hh ea ,V}h eb ]Tr [/C 1 ^, V}h 



pa r>b pc pd 
r M r M r N r N 

q 2 



(68) 



where by e&, e& we denote the length of the edges along those directions. The numera- 
tor has five Poisson brackets, whereas the denominator has the power of volume g 9 / 4 = 
t> 9//2 y 9 / 2 , which when absorbed in the five brackets, gives a contribution of v g/ 10 {h e , V} = 
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W/i e , V 1 / 10 } for each. Thus now the first term of the operator in Equation (65) has the 



10 
form: 



(io) (c em c en c ep c ea c ec y 



r> >9 q2 q2 q2 q2 

17 ^e^e^e^e/a^ 



l_ e UK e mn PTr [T'h^he^V 1 / 10 }} Tr [T J {h^V 1 / 10 } 



x Tr 



pa pb pc pd 
r M r M r N r N 



(69) 



T K h^{h ep ,V l l 10 }\ Tr [h£{h ea ,V 1 ' 1D }h eb \ Tr [h^{h ec ,V l l w }h, 

Now, the subsequent task is to find the expectation value of the operator in the coherent 
states, by first lifting the Poisson brackets to commutators. Thus now rather lengthy expres- 
sions occur, though the calculation is quite straightforward. One of the main observation in 
the evaluation of the operators which are product of quite a few terms is the fact that one 
can insert the complete set of coherent states which are resolutions of unity. Though the 
states peaked at different g are not orthogonal to each other, the overlap function vanishes 
as t — > 0, and hence the contribution to first order in t will be from the terms which are 
expectation values at the same g. In other words to quote a theorem in [8], < ip{g)\ip{g') > 
is exponentially vanishing as t — > (dfiL is the Liouville measure). 



Limit^o < ip\XY\tp > 



Limits J d/i L (g') < ^W(9') >< > 
Limit^ [< ^\ X W >< *l>\Y\if> > (70) 



+ 



< tp\X\ip' >< tp'\Y\ip > - < tp\X\ip >< ip\Y\i/j >j] 
< ip\X\ip >< x(;\Y\ijj > +O(0) 



(71) 



Thus, one now, can take the terms in (69), and then one inserts the complete set of coherent 
states so as to isolate the individual trace terms. The expression in (69), is a sum of products 
of trace terms. Thus the expectation value of the product of the operators can be broken 
into product of expectation values as follows: 

<V|yW 4 |V> 

| < ^e IJK e m ^Tr [^h^h^, V 1 / 10 }} Tr [t^^, V 1 / 10 ; 



x Tr 

C 1 
Q3! 



T K h^{h ep , V 1 / 10 }} Tr [h^{h eal V 1 / 10 }, h eb ] Tr [h^{h ec , V 1 ' 10 ^} P a M P b M P c N P^ > 



e UK e mn P < ^ Tr \ T i K i {hemj V i/ioy\ ^ >< ^| Tr \T J h-'{h^ V 1 / 10 }] W > 



x < ip\Tr 
x < ib\Tr 



T K K^{h ep ,V l ' w }\ |V >< ^\Tr [h^ih^V 1 ' 1 ^] |V > 



(72) 
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(C = 10 10 (C , em C , en C , ep C , ea C , ec ) 2 , and Q = ir ] s 2 a s 2 b s 2 c s 2 d e 2 a el Thus, once the expectation values 
of the individual trace terms are taken, a typical term in the expansion is of the form (72): 



< ^|Tr [T 1 ^ 10 - C^ 1/10 M] IV > 

~ E < ^'l Tr [( Tl Km Vl/1 ° h em)] AA IV > 



(73) 



The first term in the above vanishes due to the presence of the trace of T 1 , which is 
0. From the results of [9], one can directly replace the expression by their classical values, 
however, since the coherent states are not eigenstates of the holonomy operator, we take a 
more careful approach in our analyses of the curvature operator. One proceeds by taking 
Hab = e~ 3t / 8 [e iplTl / 2 g]AB and with K~ A b = e~ 3t / s ffie tplT I 2 \ab- The coherent states are 
eigenstates of the operators g and and hence one can extract their eigenvalues from the 
expectation value. The equation (73) gives: 



E t{ b 

ABB'CC 



i - 3t/ %B>9c'A < ^ei% ll2 V 1/W %c? J/2 \^ > (74) 



Here A, B — 0, 1 denote the SU(2) index. Using the Baker- Campbell-Hausdorff formula, 
the quantity in the brackets assumes the form: 



<V| 



p t sinhp/2 , sinhp/2 

cosh h iT Pr 

2 2 p p 



yl/W 



t sinhp/2 j sinhp/2 

cosh(j9/2) h iT Pi 

2 p p 

(75) 



where the operator p = yj P 2 + t/A. It is interesting to note that this operator never has 
the zero eigenvalue, even in the coherent state peaked at the classical area P = 0. The sum 
now reduces to 



T{ B gic9cAe- t/2 < VI 



ABC 



cosh(p/2) 



t sinhp/2 



+ e 



+ e 



-t/2 



X; TiJ^Ti^gcAe- 1 ' 2 < M 



ABB'C 



2 P 
jS inhp/2 T/1 



cosh(p/2) 



t sinhp/2 



P 



v yiot cosh P_ 

V 2 2 



2 p 
t sinhp 



P 



|V> 



sinh p/2 



P 



IV > 



~ t/2 E Ti B g BC TS cl gc'A<^(cosh P --i S ^)(P^ 

AB'C \ A A P ) 



(76) 



ABB'C A P P 

The above, clearly in the limit t — ^ takes the classical values, but which are bounded as 
f fmin- Now, the operator whose expectation value is to be evaluated is not a potentially 
divergent term, in the vicinity of the singularity, which is essentially P — > 0. Any divergence, 
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shall be in the terms including g,g\ due to their dependence on the classical holonomy. To 
isolate the potential divergent terms in the above, we simplify each of the traces in the terms 
in the above. 

The first trace term to simplify in (76) is: 

Tr(TVs) 
« ^{T I h- l e- iplTl ' 2 e- iplT '' 2 h)e m ' A 

= Tr^T'/rV^V' 74 

= Tr(T 7 e- lp/T V t/4 + e 3 '/ 4 Tr([/i, T I ]h~ 1 e~ iT ' pI ) 



= %2P 



7 sinhP^ 3t/4 _ ^ 3f/4 2q- J sincr c I J Kr Y T (r K h~ 1 c~ iTlpI ) 



P 



a 



. r sinh P 



2 .pi r ^ e 3t/4_ e 3t/4 



2a J sin u 



.UK 



a 



pK 

2i sinh P cos a 

P 



(77) 



sin a a L P M 

+2a K cosh P - 2ie LMK sin a sinh P 

a crp 

For the next few equations, we concentrate on the Trace terms one by one, which containes 

the h dependence and hence a potential divergence term. 



Tr(TVr J <7) 



,3t/4 



Tr{T I h- 1 e- iTlpI / 2 T J e- iTlpI ' 2 h)i 
Ti{hT I h- 1 e^ T,pI l 2 T J e^ Tlp, l 2 )e zt l i 
Tr(T / e" iT - p/2 T J e" iT - p/2 )e 3 * /4 - e IJK 



a J sin a r 



T E '/i" 1 e" iT - p/2 T J e" iT - p/2 



Tr(T 7 e-* T - p / 2 T J e- iT - p / 2 )e 



3t/4 



JJK 



a J sin a 



a 



cosaTT(e- iT ^ 2 T J e-^ 2 T K ) 
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- (7 L _T r (r i e" iT ' p/2 r J e" iT - p/2 T J{ ) ) e 3t/4 

The next trace term is quite similar. The other term is then equal to 

Tr(TVT J T A '#) = —Tr(T I g" f 5 IJ g) + e IJM Ti{T 1 g ] T m g) (79) 

The first term in (79) is of the form in equation (77), and the second term is as in equation 
(78). Thus the dependence on the holonomy would be precisely of the form as stated earlier. 
Clearly, all the terms of the above operator are bounded as P — > 0, which is the location of 
the singularity. The other type of the terms as obtained from equation (72) are: 

< ^\Tr [h£ [h ea , V 1 ' 10 ] h eb ] \ij> >< ^|Tr [h~} [h Cc , V l ' w ] h ed ] ^ >< ^\P a M P b M P c N P d N \^ > 

(80) 
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The terms for arbitrary a, b lead to extremely complicated terms, which we refrain from 
writing here, as the specifics are not important, and the observation of the non-divergence 
is continued for those terms also. See Appendix B for the case of the non-diagonal metric. 
What we discuss here, is where the classical metric is assumed to be the diagonal spherically 
symmetric one, and hence P^Pm = (P a ) 2 5 ab and hence, one eventually gets a simplified set 
of terms, one a factor of which would be: 



< ip\Tr 



(81) 



The first term when written in terms of QaBi would simplify to 



]Te*/ 4 <V^ 1/10 ( 



, , . . tsinh»/2 m r „ Tsinh(i»/2) 

coshQo/2 + iT J P J 

2 p p 



g BA \ij > (82) 



AB 



The first two terms simply include the trace terms of g and the last term is of the form: 



Tr(T g) = cos a cosh P 



a 



ipi 



sin a sinh P | e 3i/8 



a 



(83) 



The next term would be considerably more complicated, but the procedure would be the 
same, and the terms obtained would be of the following form: 



-Tr<^\h-y i l w hl a \^> 



(84) 



which is 



E e t/4 9l a ab < ^|e^V 1/10 [e^^]c^ ea ^[e^^]^ eaFA |^ > (85) 



A,B,C,D,E,F 

From here the trace terms would be of the form: 



Tr^s-sOe" 



9t/8 



Ti{h~ l e- %TP ' 2 e^ TP ' 2 he~ lTp l 2 h) 

Ti{e- lTP he- iTP ' 2 ) = Ti{he-^ 2TP ) 

cos aTr(e-^ 2TP ) + ^^Tr(T J r e ^/ 2TP ) 

a 



Then: 
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TrtfT'gg) = Tr(/ie- 4TP T 7 e- iTP / 2 ) 
= cos ( rTr(e- 3T - p / 2 T / )- 



a J sin a 



a 



Tl ( T J e -iTP T I e -iTP/2^ 



(87) 
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Next: 



TrtfgT'g) = Tr( e - 3iTP/2 /iT') 

= Tr(/iT 7 e- i3TP / 2 ) = cosaTr(T 7 e- j3TP/2 ) + ^^Tr(T J TV* 3TP/2 ) (88) 



Finally: 

Ti^T'gT'g) = Ti(hT J e^T 1 e~ tTp l 2 ) 

= cosaTr(T J e- lTP T I e- iTP / 2 ) + 

The second term in Equation (39) is of the form: 



a K sin a 



a 



^ T K T J e- iTP T I e- itP l' 2 ) (89) 



-yfq \K ab K cd K ef K gh q ae qf c q b °q° h 



(90) 



This has to be written in terms of the holonomy and the momentum, using the expressions 
in (50), which is: 
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^■b&dfi j '&hS e a S 'e^S e c S ed S ee S e ^ S eg S eh 



Trih;^, V]h eb }Tz{h^[h ec , V]h ed }Tr{h- e t[h ee , V]h ef } 



x Tr{/C> e9 ,V> e9 } 



no r>b jdc r>d r>e r>f r>9 r>h~ 
r l r I r J r J r K r K L r L 



(detef 



a\2 



Cf)C(iC/Ch,S ea S e() S ec S eij S ee S e ^ s eg s eh q^^ 



Tr{h; a \ [h ea , V]h eb }Tr{h- e ^[h eo , V]h ed }Ti{h^[h e ^ V]h ef } 

(91) 
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X 
X 



_ C a C c C e C g C m C n ,C n C n ,C p C,v° 1 eIJKe Tl , n ,V^°]}Tr[T J h^[h en , V^]] 
e b e d efe h s a s b s c s d s e SfSgS h (3!) 2 

Tr{T K h^{h ep , V^V^'e^'Tr^'h^ [h^, V^}}Tr{T J ' h^{h en „V^}} 
Tr[T^> v , ^ 1/10 ]] [Tr{K a \ [h Ba , V^}h eb }^{K}[h^ V^^iK^, V^}h ej } 
Tr{^K 9 , V^°)h eh } P^jPjP^PlP, 



(92) 



The evaluation of the expectation value of the above operator will follow in the precise 
way as discussed above. The expressions would correspond to the evaluation of (81). Thus 
now all the terms which appear in the expectation value of the curvature operator have 
been expressed as functions of regularised variables which are completely bounded. The 
curvature operator is a function of a particular vertex and the edges meeting at that point. 
If the vertex is limitingly taken to r — > or P — > in the regularised variables, then the 
curvature is always finite if the edge length is kept as non-zero. All the expectation value 
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of the operators are finite. From a analysis of the constants in front of the expression for 
Equation (69) and in Equation (92), the terms go as tl0 l/ 2 3 , where e is a typical edge length. 
Even if all the edge lengths are infinitesimal, the divergence reappears iff t — > 0. The reason 
is that, the expectation value of the volume operator, V 1 ^ 10 contributes from each term, to 
give finally the expression po^g~/g - Thus taking the edgelengths to also takes the Volume 
to zero in the numerator. The divergence is prevented absolutely by the observation that 
there is a minimum area as long as there is a non-zero t. This resolution of the curvature 
singularity is a calculation of the expectation value of the actual operator in coherent states 
and different than calculations of the inverse scale factor in cosmological theories [16]. 
To summarise: 

1) The Kretschmann scalar ^fgR 2 is taken and it's expectation value evaluated in the given 
coherent states. 

2) The operator is written in terms of the extrinsic curvature (65), as by construction, the 
coherent states are peaked on a spatial slice whose intrinsic curvature is 0. 

3) The terms in the Extrinsic curvature are regularised in terms of the holonomy and the 
dual momenta as in equation(50). 

4) The calculation of the expectation value of the operators gives finite terms and a propor- 
tionality to powers of 1/t. This shows that the singularity reappears as t — > 0. 

V. THE APPARENT HORIZON 

We now proceed to the question of origin of black hole entropy. Clearly we are dealing with 
a single spatial slice here, and have not tried to follow the evolution of the black hole space- 
time dynamically. The coherent state has also been constructed to obtain the geometry of 
a spatial slice. So, one cannot talk about global quantities like the Event Horizon, and even 
'Isolated Horizons' and try to obtain boundary conditions on the coherent state wavefunction 
by studying their pull backs on the given spatial slice. Hence, the relevant quantity to study 
is the Apparent Horizon, where the knowledge about the intrinsic metric and the extrinsic 
curvature of the slice are enough to determine the existence of a trapped surface of the 
equation. Thus, the geometry of the slice is built using a coherent state, such that it includes 
the apparent horizon, and then one can proceed to integrate out the wavefunction which is 
inside the apparent horizon, and obtain a suitable 'entropy of the apparent horizon'. In the 
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following discussion, we shall talk about apparent horizons which have S 2 topology. Given 
a coherent state, we try to re-build the classical space-time by evaluating the expectation 
values of momenta and the holonomy in these states. We recover the information about the 
extrinsic curvature and the intrinsic metric of the slices. Question is how does one know 
that there is a apparent horizon in the slice? There is a very general equation, which when 
satisfied shows the existence of a trapped surface. This is of the following form: 

V a S a + K ab S a S b - K = (93) 

Where S a is a space-like vector, normal to the 2-surface. The above can be re-written in 
terms of the variable K T a and the triads as: 

V a S a + K I a e I b S a S b - K 1 ^ = (94) 

For arbitrary value of the Immirzi parameter, this equation will show a dependence on the 
parameter. Here, as 13 = j3Ki and ^e* = (l/y/]3)e I a , ^e? — vT^e", the above equation in 
terms of the generelised variable assumes the form (the Christoffel connections which appear 
in the covariant derivative are independent of @) : 

V a S a + f>Kl p e[ - ^ *K{ "e? = (95) 

Now, we stay with the assumption that we are trying to re-build the Schwarzschild black 
hole space-time from the coherent state, and hence, there is a set of spherically symmetric 
coordinates. The apparent horizon, is a sphere or S a = (1,0,0), and one gets the following 
relation for the intrinsic metric measured, and the extrinsic metric of the same slice. (We 
resort to the classical variables or 0=1). 

-Tl - rj r + K„ - K rr q rr - K m ^ - Keeq 96 = (96) 

Now, T e '^ r , Kg g, K^rf in the classical values cancel each other at r = r g leaving the following 
equation in the radial sector, which is trivially satisfied everywhere: 

K rr (l-q rr ) = (97) 

Thus clearly imposing the apparent horizon equation on the radial edges does not introduce 
any new constraints on the radial coherent state wavefunction. One can then think of the 



31 



horizon as being formed by purely radial edges crossing the horizon, with the radial wave- 
function corresponding to 'free wavef unctions' as would be there for a spherically symmetric 
space time. 

v 




V 

o 



While trying to lift the apparent horizon equation to an operator equation, we only retain 
the terms which constrain the angular sector of the theory. A crucial assumption which 
is included in the calculation is the fact that classical spatial slice has a flat spherically 
symmetric metric, and the horizon is a spherical surface in that spatial slice. 

The apparent horizon equation now only includes the derivatives of the angular metric 
forming a difference equation for the discretised variables along the angular edges. The 
equation is: 

-It - - V% - Kief = (98) 

The Tg r = \q ee (qee),r = e^e^egj t r egj) The derivative is replaced by a difference in the 
discrete version of the same equation. The horizon, is thus now a set of radial edges, with 
a set of vertices {vi} outside the horizon, and a set of vertices {t^} inside the horizon. The 
derivatives are thus differences in the value of the metric at vertices {v\}, with their values 
which are at vertices {i^}- Thus for a representative set of vertices, the difference equation 
is obtained by subtracting the value of eej at v\ which is outside the horizon with the value 
of the inverse dreibein inside the horizon at vertex i> 2 . The inverse dreibein is also expressed 
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as a Poisson bracket first and then lifted to an operator form. Thus now 
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e r e 
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Tr 



T J hf [h e ,V] 



Tr 



T J hg l [h e , V] 
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x Tr 



T J hf[h e ,V] 



(99) 



Further, one uses the observation that /3 = F^—j3K[, where (5 is the Immirzi parameter, 
to write -fidp ^Aj K{. This enables one to rewrite the apparent horizon equation interms 
of the holonomy operator, incorporating information of the extrinsic curvature. Thus 
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The equation is then of the form : 
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1 d 
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■Tr[T 7 %] p Pf + (0 
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— Tr 



j ' ! i 



T J hf[h e ,V 1 / 2 } Tr[T J h e [h e ,V^ 2 } 



Since the first set of terms is independent of j3 no attempt has been made to write the (3 
index there. 

Fortunately at the classical level, the difference equation has a linear term in r to derive. 
In other words, the derivative is of the form — ~(r — (r + 5)) = r, and hence yields a exact 
answer for the derivative, introducing no further fuzziness in the equation for the location of 
the apparent horizon, for fine graphs. An evaluation of the equation in the coherent states 
leads to the following equation: 



E {T^Uea < M (e lTP6/2 ) CD V 1/2 (^' 2 ) DE ^ > 

ABODE 

1 d 
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< ip\Pf\ij > 
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(102) 
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Now, in the above equation, all the quantities on the LHS contain quantities which are 
functions of variables defined inside the horizon, whereas, rest of the quantities are of vari- 
ables defined outside the horizon. Now one can proceed and evaluate the traces like in the 
curvature operator case, and obtain complicated expressions on both sides of the equation. 
The solution for the equation, exists and in the v\ —>■ limit gives r = r g . Without going 
into the explicit solution, one can infer some qualitative features which are sufficient for the 
discussion of the entropy. 

1) The quantisation introduces a 'fuzzyness' with the location of the horizon, which is pro- 
portional to the semiclassicality parameter. 

2) the classical values of gg, for the edges meeting at the vertex v\ are correlated with the 
same for edges meeting at the vertex v-i- In other words, 

9e(v 2 ) = 9e(v2){ge(v-o) (103) 

3) The apparent horizon equation is completely independent of the holonomy along the radial 
edge which crosses the horizon. 

4) The Volume operator at the vertex t> 2 contains information about the horizon Area P#, 
but we consider that as an independent variable, which can be changed, and the rest of the 
variables will change accordingly. This will imply a sum over all possible graphs interpreta- 
tion for the entropy, which we also try to discuss at the end of this article. For a fixed Ph, 
the above apparent horizon equation, gives the graph degrees of freedom inside the horizon 
as a function of those outside the horizon. These will reflect in a correlation between the 
wavefunctions inside and outside the horizon. To emphasise again, the radial edges which 
cross the horizon remain completely free, with no restriction on them from the apparent 
horizon equation. 

Clearly this results in a set of wavefunctions for the coherent state peaked at the edges 
which are correlated with each other. Now, one has to ensure that this results in a entropy 
proportional to the area of the horizon. We show in the final density matrix calculation, it 
is the radial edges which constitute the 'entropy' of the system by counting the number of 
ways to build the apparent horizon 2-surface area. 

In more general terms, very crucial assumptions which went into the determination of 
Equation (98) to have a separation of the angular and radial components was the assumption 
that there are no cross terms in the extrinsic curvature, or K r eu\ = even at the operator 
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level. This can be ensured in terms of the gauge connection by the condition 

K ab = (Ti - Ai)e n = 0, fora ^ b (104) 

and 

q ab = E a E b = ofora ^ b ( 105 ^ 

(Which is the diagonalisation condition of the metric). 

However, these assumptions are not drastic, as any deviations from the above due to 
quantum fluctuations will be proportional to t, and hence ignored in the correlators being 
considered here. Thus, the apparent horizon is located by measurements in the angular 
coherent states, and the radial coherent states associated with the radial edges, induce 
the apparent horizon with area. Now, we shall address the question why the degeneracy 
associated with the induction of the Area of the horizon can be called the entropy of the 
black hole. Hence, in the next section, we first review, how entropy arises from a density 
matrix, and then derive the density matrix associated with the coherent state wavefunction. 

VI. THE DENSITY MATRIX CALCULATION 

The density matrix by definition for a given configuration in a state \ip > is of the following 
form: 

P =\iP><iP\ (106) 
The expectation value of any operator can be evaluated in this Matrix as : 

< A >= Tr(Ap) = Tr(A\ij >< $\) =< V>LA|<0 > (107) 

The usefulness of the density matrix is particularly evident, when the system is in the 
product of two Hilbert spaces, and one defines a 'reduced' density matrix by tracing over 
one of them. As a illustration, let the system be in the state: 

\ip >=J2 d im\i> \m> (108) 

im 

where the orthonormal states \i > \m > have a product structure, with \i > and \m > 
belonging to two different Hilbert spaces. One can define a reduced density matrix, by 
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tracing over the \m > states. 



Pred Tv r , 



d im \i > \m > d* n < j\ < n\ 



= YjY^d im <r jn \i >< j\5 mn 

im jn 

= J2 d imd* jm \i ><j\ (109) 

ijm 

Interestingly, the diagonal elements of the density matrix are the probability of finding the 
system in a state \i >, given any state in the Hilbert space which has been traced over, or, 

Predii = 2^ d im d* m (HO) 
m 

and hence the condition which naturally arises is: 

Tr(p rcd ) = 1 (111) 

If there are no correlation between the Hilbert spaces, in particular if the coefficients factorise 
as di m — did m , then the resultant reduced system is still a pure system. This can be checked 
by the following condition: 

p 2 = p the State is Pure, p 2 < p the state is mixed (H2) 

As reported in the previous section, the Coherent state for a black hole has precisely the 
same structure of product over Hilbert spaces. Now, is there a observer for whom only a 
part of the system is measurable? At the classical level, (one believes that the measuring 
instruments are classical), for an outside observer, the inside of a black hole does not exist. 
He/She does not receive any signal from inside the horizon. For such an observer, there has 
to be a reduced density matrix. However, as the coherent state is written in the covariant 
form, there appear no information of the continuity conditions which must be , imposed at 
the vertices. Thus, the classical labels of different edges are not really independent, but must 
satisfy a continuity condition at a shared vertex. This introduces correlations even at the 
classical level. These continuity conditions must be imposed from Einstein's equation. Here, 
we discuss only the classical correlation introduced due to the apparent horizon equation. 
Thus even at the classical level, and one cannot integrate over the edges inside the horizon 
independently of the edges outside. There is also a method to quantify the relation. 
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A. Correlations 

Classically, we have defined the graph degrees of freedom for the black hole, as the 
holonomy and momentum along one-dimensional edges and corresponding dual surfaces. 
Though, we have talked about recovering the 'classical action' principle for these variables 
on the individual edges, we have yet to examine the continuity conditions at the vertices, 
since the classical metric itself is continuos throughout the manifold. These conditions will 
ensure the existence of correlations for the states defined on the entire manifold. At the level 
of the gauge invariant 'coherent state' wavefunction, the Gauss constraints are satisfied with 
the introduction of the Clebsch Gordon coefficients at the vertices. There is a discussion 
on the Gauge invariant states in [8], and the peak behaviour of these states. However, we 
are mainly interested in quantifying the correlation at the vertices, and in the semiclassical 
limit. This limit is difficult to evaluate from the exact form of the coefficients. We shall 
discuss this eventually, but for the time being, we discuss the 'macroscopic correlations' 
which are manifestations of the continuity of the classical metric, and at the horizon, the 
apparent horizon equation. 

For the classical observer, who makes measurements in the Schwarzschild coordinates 
t, r, 6, 0, while confined to the slice r = c slice, the spatial slice ceases to exist after the 
horizon. Though the transformed discrete variables will be the same set of variables of 
the observer, as 4^ = 1, but as he/she approaches the apparent horizon, for her/him, 
the coordinate time measured in t will be running to infinity. (This is the same graph, 
same slicing, but measured in the coordinates of the outside observer, and not the r, R, 8, 
frame.) For the outside observer, the coherent state wavefunction inside the horizon is 
irrelevant information. He/she is forced to deal with a reduced density matrix instead. The 
expectation values of all operators for the outside observer pertain to the density matrix 
now. What does tracing over the 'Inside edges' mean? The gauge-invariant wavefunctions 
are correlated at the vertices by the Clebsch- Gordon coefficients. And it is a difficult task 
to evaluate the exact nature of this correlation, and how the coefficients influence the peak 
of the entire coherent state. However, the fact that the elements g themselves are classical 
solutions restrict the wavefunctions. As mentioned in the previous section, the apparent 
horizon equation clearly correlates the classical elements across the horizon. We now try to 
quantify the correlation across the horizon. Let us take the set of radial edges which cross 
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the horizon. These end at a set of vertices, and these vertices are the starting point of a set 
of edges, now inside the horizon. If the radial edges which cross the horizon induce the set 
of classical elements g = e 4P ' Ti//2 /i e , then the areas induced by the edges inside the horizon 
also get determined from the apparent horizon equation. 

The apparently uncorrelated wavefunctions are now naturally correlated. Let us take the 
specific coherent states which are peaked over the edges sharing a vertex, one inside the 
horizon, and one outside the horizon. For simplicity, we will be suppressing the azimuthal 
direction. The edges are e ro , tQ o eQi o ,eu-,^r i -,^e i -,^e' i - Where the first three share a vertex V\ 
with the radial edge e# crossing the horizon, and the last three sharing vertex v% which 
is inside the horizon and end point of e#. To avoid overcounting of the edges sharing the 
vertex v\ and v 2 , one takes the set of angular edges which are ingoing at a given vertex. The 
coherent state is thus a product of the coherent states for each of the above edges. 

II II j' s 



x e 



-t^fsU's+l)^ > % > > b . > {h > (n3) 



The j's notation denotes all the spin eigenvalues j a , ....ki. 

Now from the apparent horizon equation, classical values of gg i = ge t {9do)i an d hence 
the the wavefunctions cannot be integrated as independent, for the different edges. The 
complete wavefunction should be of the form: 

\^(v 1 ,v 2 )>= ^Hoyimmlho} > \3h > \j{i} > (H4) 

j{0}JHj{I} 

The index for the horizon edge is written separately as they occupy a special status, as they 
are more or less free wavefunctions which are factorisable from the rest of the wavefunction. 
j{o} denotes all the spins labelling the outside edges and jyy denotes all the spin labels 
for the inside edges. The horizon wavefunction however is determined by the area Pjj 
for the horizon bit that the edge induces. Infact, it is precisely the Pjj of the horizon 
edge which can be considered as a variable and adds to the entropy of the horizon. The 
form of the coherent state wavefunction remains the same as above, except that now the 
labels g V2 of the wavefunctions inside the horizon, are functions of those outside. The exact 
functional form will depend on the solution from the apparent horizon equation. Moreover, 
even so, the correlation is difficult to quantify in a entropy calculation. To enable the entropy 
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calculation to first order in the semi- classical parameter, we encode the correlation in a 
conditional probability function f(g ra , g$ \gm 9e^) — N, (N is a constant), when the values 
are as determined by the apparent horizon equation or is otherwise. Given the formalism 
discussed in the paper, this is the only way one can introduce a correlation in the wave 
function. The 'classical correlation' is encoded in the expectation value of the operators, 
and hence the labels g e . The functional form of the wavefunction and it's dependence on g e 
remains unchanged, as obtained from the coherent state transform. 




|*(vi,ua)>= E i ) j {0} j H j {1} (nJ)\ho}>\3H>\j{i}> (H5) 

i{0}3H3{l} 

Here, the coefficient tpj {0} j H j {I} now includes the correlation function f. Now clearly the 
above coherent state is written for the edges which are connected across the horizon by one 
horizon edge. In otherwords they share the vertex v\ and v 2 which are the initial point and 
final point of the radial edge crossing the horizon. Thus the density matrix constructed of 
the above state can now be written in the form: 

\y(v 1 ,v 2 )><y(v 1 ,v 2 )\ (116) 

The full density matrix for the entire horizon would comprise of 

p = J] |*(ui,u 2 ) >< *(«i,v 2 )| (117) 

V1,V2 

Where the product is over all the vertices comprising the immediate outside and inside of the 
horizon. The tracing is now done over the edges which share v 2 , or the internal edges. The 
reduced density matrix then has the following components (following the same procedure as 
in Equation(109)): 

Pvi j{0}3H,j' H j {0 '} = ^2$3{0}jBj{I}' l i , j{0'}iHHl} 
{1} 
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Explicitly in terms of the coherent state wavefunction, the sum results in terms of the form: 

n . . ., . _ y ^{Q} *3h n f H 71 Ho'} Km (fl) I 2 1 f . \2 p -tj {1} (j {I} +(LU q \ 
Pvt Ho}3 h,J hH o' } ■ - ||^||l/2||^||l/2||^||l/2||^||l/2 II^U \JH0Wi}\ C ^ LJ > 

x e -*j{o}0'{o}+i)/2 e -tjH(jH+i)/2 e -ti^(j^+i)/2 e -*i{ 0} (3'{o}+ 1 )/ 2 (120) 

The result of the summation of the indices jyy is a Gaussian, peaked at the classical value 
of gj. (Note j c \ for any of the edges outside, inside, or at the horizon indicates a peak value 
which equates to classical parameters.) 

The density matrix which is reduced now has the following form 

_ 7r J{Q} 7r J{Q} -tj {0} (j {0} +l)/2 p -^ 0} (i{ 0} +l)/2 ^H^g t j H (j H +l)lf . 12 
Pvi 3{0}3Hf {QI} 1 1-01 I IMI l-0{0}J{Icl}l 

(121) 

Due to the | j-0 1 1 in the denominator which contributes with a factor proportional to 
e~ p2 /*e _ */ 4 sinh(p)/p, the 

(Note cl for any of suffixes for edge labels outside, inside, or at the horizon indicates a 
peak value which equates to classical parameters.) 

offdiagonal elements are all damped exponentially in the semi-classical regime where 
t — > 0, and P is very large. For the diagonal elements this factor is assimilated into a 
Gaussian which becomes a delta function in the semiclassical limit. Infact as determined in 
[8], this implies 

7f-7T- / n j (m/j- i? P 3 /P) 2 j (n/j- L P 3 /P) 2 r f , +1/2U _ Pl 2 

W e -^0+i) „ e 2 (i^p 3 /V (i-^^/n 2 e - [b+1 T 1 (122) 



The interesting aspect of this is that when the semiclassical parameter t — > 0, the functions 
in (122) tend to delta functions, each peaked at appropriate values of the classical variables. 
Thus, the 'density matrix' reduces to a set of delta functions: 

= S(j{o}t, P {0} )S(m {0} t, P( 0} )S(j H t, Pn^imHtP^lf^J 2 (123) 

The density matrix has non-zero values for only a set of matrix elements, as determined by 
the delta function equations. The non-zero elements are Pv 1 j {0)cl j Hcl m Hcl ,j Hcl m Hcl j {0}cl , however, 

the rriHci = —]Hc\ Jflci- At the horizon only Ph is fixed, and hence the entire range of m# 

is allowed. This is the origin of the degeneracy of the density matrix, and the entropy of the 
entire space-time. After all the discussion on the coherent states and the density matrices, 
the entropy is finally the number of ways to build the horizon area. Inclusion of all the 
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edges comprising the entire graph in the initial density matrix would not contribute to the 
entropy, as they would belong to tensor product of Hilbert spaces, and give unit norms in the 
trace. However, one must emphasise that this is a 'semi-classical' result, valid at 0th order 
in the semi-classical parameter t. If one includes the higher order corrections, one would get 
a 'quantum entropy' and the importance of the above result for an explicit expression for 
the density matrix is in the obtaining of the corrections to black hole entropy. From (123), 
and Tr(p) = 1, it follows that |/j 0clJIcl | 2 = 2 j H \+i ■ This apparently 'trivial' factor / still has 
non-trivial information about the classical correlations, as now jo c \ and jj c \ are determined 
by the classical labels of the edges at the vertex v\ and v<i- Hence, if the jo c \ and ji c \ were 
anything other than those determined by the apparent horizon equation, then / would be 
zero. Now the full density matrix will be a product of all the density matrices of the vertices 
immediately outside the horizon. Thus the Entropy would be: 

Sbh = -Tr[plnp] = -Tr[]]_p v lnfljp,,)] = X>(2j„ + 1) (124) 

V V V 

The additional constraint is 



p 



In case all the area bits are set as equal j v = j s , one obtains, 



N(J S + \) = ^f (126) 

and Entropy 

SB - = u^m HV - +1) (127) 

This obviously differs from the Bekenstein Hawking entropy due to the dependence on 
Spin. We will subsequently discuss the Immirzi parameter which will be adjusted to give 
the appropriate entropy. 



B. Sum over all possible graphs 

In the previous discussion, the entropy is shown to arise due to the degeneracy associated 
with the number of ways of inducing horizon area, given a graph. The graph is such that only 
radial edges cross the given horizon. Now, this assumption can be somewhat generalised 
to include a sum over all possible similar graphs, which would imply a difference in the 
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number of edges crossing the horizon. This then would include a sum over all possible areas 
induced by a horizon edge, though the entire formalism of the previous section would remain 
absolutely same. The coherent state formalism is useful is studying this sum over all possible 
Ph for the reason, that, there exists a Liouville measure in the classical phase space, with 
respect to which the states are overcomplete. Now, when the density matrix for a particular 
set of vertices is defined as 

p vlV2 = |*(vi,u 2 ) >< V(vi,v 2 )\ (128) 

The apparent horizon equation fixes the value of the correlated wavefunction ipj {0} j H j {I} at 
a particular value, as all the classical labels are determined. However, if now one chooses 
to integrate over all possible areas induced at the horizon, then, the classical labels for 
the ge a and gg i are also going to change accordingly. To quantify the above, one can thus 
write the conditional probability as now a function of Ph, with a j (ge , ge^Pii) = N, if the 
classical labels are as obtained from the apparent horizon equation, or zero otherwise. Once 
the density matrix is obtained, then the diagonal terms give the probability of finding the 
system in that state, given any internal state. Now, this probability function is multiplied 
by the 'overlap function' or the probability p t {gn, 9h)i with respect to the Liouville measure. 
The probability function gives the probability of finding the system at the phase space point 
g', when it is in a state ^(gn)- 

p\g H , g' H )dQ = Ag H , g^e^e^ 2 ^ /2 ^ H d 3 P H d^ H (129) 

Where d\xn is the Haar measure. Now, finally, the probability and the correlation function 
are independent of the holonomy, and so, the Haar measure integrates out to 1. Also, in 
the t — > limit, the ovelap function is almost 1. Thus, one retaines now only the Liouville 
measure in the variable Ph- Now the measure d 3 Pn can be converted to P^dPn sin Qd<dd§, 
where the angles G and $ are in the internal 3 space. Integration over the angles would 
achieve a averaging over the components of the momentum, with it's length as fixed. 

The diagonal terms in the density matrix dependent on the horizon wavefunction are of 
the form: 

^) mm e-^) e -^/* e -'/*^? sin QP 2 H dP H (130) 
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The expression for 7Tj(H 2 ) mm is then the following: 

I tU 2 \ 12 V- (j + m)\{j - m)\ r 2 2 2 p 

^(tf U| = E (j _ m _ /)!(j+m + , (/!)2 [cosh P„-cos 6smh P H 



cosh P# + cos G sinh Ph 
cosh Pff — cos G sinh Ph 



sin 2 G sinh 2 Ph 



cosh Ph — cos 2 Q sinh Ph 



(131) 



Using this back in the equation, the above leads to a integrand which can be reduced to: 

E E E M H™ cosh 2 ^' P H [tanh p Jf ]2(7-*-*72) ^y-m-l-k^-m-l^m [cog Q] 2(j-m-Z-/=) gin2; cog2m - 
l k k' 

(132) 

with 

M (j + m)\(j - m)\ 

Jlm ~ (j-m- + m + /)!(/!) 2 1 J 

The integral over Q is /J r sin 2Z+1 6 cos 2 j~ 2l ~ 2k ~ k ' QdQ and hence yields a beta function B(2l + 
2, 2j — 21 — 2k — k' + 1) , with the additional restriction that A;' = 2n, such that the sum is 
non-zero, the sum now reduces to 

E E E M Jim cosh 2(j) P H [tanhP H ] 2 °'~ fc ~ n) {-l) j ~ m ^ k C{- m - l C^ B{2l+2, 2j-2l-2k-2n+l) 

l k n 

(134) 

with 

M (j+m)!(j-m)! 

M * m (j-m- + m + /)!(/!) 2 1 j 

Defining k + n = q one gets: 

EEE^ m cosh^P, [tanhPH] 2 o-) ( -iy-m-«- fc ^ ca £(g + 2)r(2j - 21 - 2 g + 1) 

(136) 

This is thus the expression for the density matrix element, and a integral over the Liouville 
measure will give the 'sum over all possible graphs' entropy. 

There is a interesting observation due to the e~ p H jt term in the Liouville measure, this 
is discussed here. The calculation simplifies due to the observation, that classical Ph = 
sin 6 6<f), if one takes a = r 2 , the only possible area scale in the system, and hence independent 
of the black hole horizon radius! The value of Ph should vary from to oo, depending on 
the number of edges crossing the horizon (increase in number of edges decreases the value 
of Ph)- Also, due to the presence of the e~ p H^ term in the measure, the maximum value of 
the density matrix elements are concentrated around Ph = 0. In the regime where Ph ~ 0, 
the / = term in the sum dominates, and the sum trivially reduces to 1. This concentration 
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of the possible Horizon area elements near a (j c i + l/2)£ ~ value is in agreement with 
previous derivation of black hole entropy as arising from j = 1/2 spin elements, albeit in 
different formalisms [19]. The density matrix elements then are of the form: 

Pii = \f\ 2 (137) 

Setting J2 m \ f\ 2 — 1) one obtaines \ f\ 2 = l/(2j c / + 1), and hence the correct value of entropy 
is recovered. This limit is in the situation where the number of edges crossing the horizon 
are large, and quite opposite to the regime discussed in the previous fixed graph calculation. 
However, in both the regimes the same entropy law, proportional to area of the horizon is 
recovered. Further, in this regime, the offdiagonal elements of the density matrix will also 
be non-zero, and provide corrections terms to the entropy. 



VII. ENTROPY AND IMMIRZI PARAMETER 



This factor prevents the above entropy counting from being the exact Bekenstein-Hawking 
entropy. However, what has been excluded from the above is the fact that the area spectrum 
is {jci + 1 /2)/3t, for different quantisation sectors of the theory. Since the constraint is 
N(j s + %)lp = Ah, where Ah is the classical area of horizon, which remains unchanged for 
the Immirzi parameter related variables. 

N = MTW) H2j - + 1) (138) 

Thus it seems that if /3 = wij+f) wom d be the correct choice for the Immirzi parameter. 
However, this is rather dependent on the choice of graph, and the area it induces on the 
horizon. In the most general situation, this would include a sum over all possible graphs, 
and the coefficient of the counting would equate the value of the Immirzi parameter. Here 
the counting would be different from that given in [18] and hence the value of the Immirzi 
parameter would be considerably different, from the prediction there. 



VIII. ENTROPY OF SPACE-TIMES WITH APPARENT HORIZONS 

Here we summarise the results in the previous sections for the entropy of the Schwarzschild 
black hole, and show how they can be easily extended to include spherically symmetric 
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apparent horizons. The key assumption of the entire article has been that the classical 
spatial slice, at which the solution is peaked has spherical symmetry. In otherwords, the 
metric on the spatial slice has the form 

ds 2 = dr 2 + r 2 dVL (139) 

This also presumes that the intrinsic metric is flat. The information about the curvature 
of the entire slice is contained in K a b or the extrinsic curvature for the entire slice. The 
presence of the apparent horizon equation is obtained as a solution of the equation 

-Tl - r*. + Kiel + Jtje} = (140) 

Where Kgg, are arbitrary functions of the radial coordinate. As before the graph is 
taken along the radial, and angular coordinate lines, forming discrete lines. The classical 
space is thus sampled in terms of discrete variables h e , P/, which are then combined to form 
the classical label on which the coherent states are peaked. Since the apparent horizon 
equation places no conditions on the radial edges, a similar construction is done, where 
horizon areas are induced only by the radial edges. The apparent horizon equation then is a 
difference equation linking the vertices which the radial edges connect. The coherent states 
are then correlated across the vertices. The density matrix and the subsequent derivation 
of entropy thus follows automatically. This would imply e.g., any spherically symmetric 
apparent horizon will have entropy proportional to the area of the horizon, this would then 
automatically include cosmological horizons, also. The addition of a cosmological term in 
the lagrangian does not change the definition of the kinematical variables of the theory, and 
hence all the discussions in the previous sections of the coherent states does not change the 
derivation of entropy. 

IX. CONCLUSION 

In this paper, the main aspects of black hole physics emphasised are: 

1) Black holes are classical solutions of Einstein's equations and hence would correspond to 
a semi-classical state like the coherent state in a quantum theory. 

2) In a exact quantisation of the theory, it is difficult to identify a complete 'black hole' 
state, though, states with trapped surfaces in reduced phase space quantisations have been 
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obtained in [17], and also the isolated horizon boundary condition basically identifies the 
horizon by imposing the boundary condition on the physical states of the theory. 

3) The coherent states defined in [7] on the other hand provide a complete formalism for the 
state which would correspond to a given black hole solution. 

4) The given coherent state which is peaked in both the momentum and configuration space 
variable is also a state where both the 'area' (intrinsic metric) and the extrinsic curvature 
(holonomy) can be measured with minimum uncertainty. 

5) Thus the measurements in the coherent state allows one to build a entire spatial slice of 
the given space-time, and hence locate the apparent horizon and the singularity uniquely. 

6) The information about the space-time is however encoded in regularised variables of the 
holonomy defined along one dimensional edges of finite length, and corresponding momen- 
tum induced on dual surfaces of the edges. This regularisation in case of a singularity at 
one point effectively regularises the same. The holonomy which has information about the 
curvature never diverges and takes oscillatory values from -1..1. 

7) A clear uncertainty principle prevents measurement of area zero, and hence, the singular- 
ity in the curvature operator is resolved mainly due to the existence of a minimal area. 

8) The Coherent state is a tensor product of coherent states over edges, and hence one can 
trace over the edges inside the horizon. 

9) To encode the correlation of the classical variables across the horizon, one finds that the 
apparent horizon equation is a function of the variables outside the horizon, as well as those 
inside the horizon. 

10) The wavefunctions which are functions of the classical labels are thus correlated, and 
lead to a entropy of the resultant density matrix obtained after tracing the wavefunction 
inside the horizon. It shall be very interesting to calculate the corrections to entropy, and 
whether they have a logarithmic behaviour as calculated in [22]. 

There are quite a few questions still to be completely understood. 

1) What exactly happens for a general spatial slicing, i.e. those without the flatness as- 
sumption? 

2) Are there more coherent states than the ones stated here, and what would the expectation 
values of operators in those states be? 

For immediate future problems: 

l)The above formalism can be generalised to arbitrary black holes, those with charge and 
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angular momentum. 

2) The regularised curvature operator in terms of holonomy and momenta can be measured 
in other states, not necessarily semi classical. 

3) The semiclassical derivations of black hole thermodynamics here can be used to obtain a 
approximate temperature using the density matrix, and a evolution in time. This is work 
in progress. 

To summarise, the semi-classical nature of black holes have been discussed. The quanti- 
sation, is non-perturbative, and includes the full degrees of freedom for gravity. It will be 
interesting to verify the Hamiltonian constraint completely for the coherent states, right now 
the action of the Hamiltonian on these states leads to states of very small norm. There is a 
discretisation of the Schwarzschild space-time using dynamical triangulation techniques in 
[20], it shall be interesting to obtain the entropy in that formalism. 
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Appendix A: The Quantum Action Principle 

The Quantum action principle is postulated as 



Applying it to the restricted set of Coherent states, one gets the following expression: 




(141) 




(142) 



(dfi(h) corresponds to the Haar measure) Now, to quote [8], the state 



(143) 
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As calculated, 



Xj(gh ) 



(144) 



A -A- 1 



with 



A = x + Vx 2 — 1 



(145) 
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where x = cosh(P/2) cos 9 + i sinh(P/2) sin(6 l ). Thus 

^foA ') = -j [ Xi (9h + ft + 1) A — A -1 J (146 » 

Putting this back into the scalar product, one obtaines: 



d r A 

< vi^-i^ > = y dn(h)- 

1 d 



(147) 



The first term in the above integrand clearly is proportional to the probability density, which 
as t — > 0, is peaked at the value of 9 = 0. We investigate the second term there, to see if 
the sum can be converted to a convergent expression. The technique used there are as in 
[8]. We concentrate on the sum first: 

\ (2j+l) , \-(2j+l) 

E(2j + 1) v^ + i)/2:L_±A_ — 



^2 e -(t/8)(2i(2j+l)) 



A-A 



e (2j+l)lnA _ e -(2j+l)lnA 



1 



T £ 



2 -t/8[n 2 -8nlnA/t] 



A-A 



n e 



1 oo 

c */8 1 V" n 2 e -[(ns) 2 -2(ns)z]/2 

A — A 



n=— oo 



with n = 2j + 1 and s 2 = t/4, z = 

This is converted by the Poisson Summation formula to a convergent series by using the 
following property 

E /M = — E /( — ) ( 149 ) 

n=— oo n=— oo \ * / 

where / is the fourier transform of f(ns). 
Here, 

S 2 f(ns) = {n S f e -\{ns?-lns) Z ]/2 

The fourier transform of which can be easily evaluated as the integrals are simple Gaus- 
sians. 
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f(k) 



2 + (z- ikf 



e 2 



Uz-ik) 2 



Substituting the same in the Poisson summation formula, one obtaines 



(151) 



n=—oo 



2. - 



2 + [z-i- 



.2tt 



As the sum is computed, in the t — > limit, only the following terms are relevant 



^7T 



1 [cosh _1 (x)] 2 J 



A-A- 



^[oMh-^x)] 3 ^ 



[cosh (a;) 



(152) 



(153) 



with In A = cosh (a;). Substituting in the expression for the action, this term in the action, 
one obtaines the following set of terms: 



' resquantum 



dfi(h) 



(A + A^ 1 ) cosh _1 (x) 



A-A- 



p\h) 



where p t (h) = ipip/ 1 | 2 . Now, 



.r 



A 



A A- A 



-i 



(154) 



(155) 



with 



x 



sinh 



cos 9 + sin 99 cosh 



l>>9> + Pitf , 
+ i — sinh 



P9 



sin 9 — i 



P 
~2 

pjQj 

P 2 9 2 



P9 + 9P) sinh 



sin 9 



+ z- 



pjgj 



P 



cosh — P sin 9 + %■ 



piQi 

w 



p 



sinh — cos 99 



(156) 



Substituting the above in equation (147), one obtains the terms which survive after the 
integration of inner product, (one replaces the terms with 9, 9^ with 0, as the coherent state 
is peaked precisely at these values), 



also, 



ijj d 1 
-^ = 1 



E(2J + 1) 



P— + i9 j P j 
2 



-tj{j+i) 



UH 2 



one obtains: 

i 



d 

dr 



^ 3 



(2j 



'£2?+! + £-(2j+l)\ £ + £ 



-1 



(157) 



(158) 



(159) 



19 



Given the fact that £ = coshP + sinhP, the above can be calculated, and in the P — > oo 
limit, it gives: 

d n./.ii2. t ^(2j + 1) 2 f e 23+1 + r ( ^ +1) l e - tj( j+i) (160) 



The second term in the above yields in the P ^> 1 limit 



t 3 / 2 sinh p , 



OO 



E 



e 



(z-ik) 2 /t 



-oo 



(2 — i/c) 



t 3 / 2 \ t 

Thus finally, in the limit P — > 00 









IMI 2 = ip 




-t) 









161) 



--PcothP » --P (162) 
The first term needs a Poisson resummation formula implemented (z=P): 

— ^— £(2j + l) 2 + r W+1) }^' +1) (163) 



(164) 



e ^/t 1 + — (165) 



1 d ' ^ l P (166) 



dr n " 2 VP t/ 2 
Interestingly, substituting the above in equation(147), one obtaines the cancellation of 
the first term in Equation(157), and finally one obtaines the expression for the classical 
action. 

Appendix B: Trace terms in the Curvature Operator 



Since the trace terms take rather long expressions for certain cases, they are described 
in the Appendix here, though the physical interpretation for them is exactly same for each 
of the terms, and the conclusion is that the expectation value of the curvature operator is 
bounded. In the evaluation of the terms which are of the form 

< ^\Tr[h^[h ea ,V]hM > (167) 

The following traces occur: 

pb papb pb 

Ti(glg a g b )e- 9t/8 = 2 cosa 6 coshP a cosh — + 2 cosa 6 -^ sinh P a sinh — 
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coshP a -2«^P; 



j jsino-fcSinhP 6 / 2 

" 2 ^ P ^ P^~ 

r\ ijl P a pbL sin cr 6 sinh P b /2 sinh P a 



a sin a b sinh P a /2 



(T 6 P a 



cosh(P 6 /2) 



a b P a P b 



2iP?a 



I u a 



a ^ jS i nhP a sma a ^ IJK 
P a a a 

sma a KMN M B sin a a sinh P b /2 
2a^ cos cr fe cosh(P /2) 2ze er™ "n cos a f> Tn, 

® a a* 

t , h , . sin (Jh wmm m m sin a?, sinh P 6 /2 

-lal cos a a cosh(P b /2) - 2ie K MN a ^p^ cos a a 



+ e &b a a cosn(r^ /zj 



6 sin cr 6 sin a a 



'b 1 b 



(T b P b 



Ob<Ja 



(168) 



The other trace term is 

Tr [gl Tl 9a9b 
which can be simplified to 



e _ 9t/8 = ^ ^ lPa T/2 K) ( e ^P b T/2 hb) { j l -l e -iP a T/2 )T I-\ (16g) 



Tr 
= Tr 
-2i Tr 



e-^Haih-^e-^he-'^T 1 ] + Tr \( e - lP ^ 2 h a [g b , h^^T 1 



hbe -iP aT /2 T I e -iP a T/2 e -iP b T/2-\ + 2Tr ^-^2^-^/2 



ijk i jsma a sma b 

e 

<?a &b 



h n T k h b e- w « T ' 2 



cos cr&Tr 



iP a T/2rpI e -iP a T/2 e -iP b T/2 



+ 



0V, 



2e 

<?a Vb 



T J e -iP a T/2 T I e -iP a T/2 e -iP b T/2 

Tr{cosa a [ e -^/2 r * e -iP T/2 



erf sin a 



(T 



T L e -iP b T/2 T K e -iP a T/2 



} 



2i 



cos a a cos CT(, — cr a a"6- 



sm a a sin (75 



0"a 0"6 



Tr[T K e- iP " T / 2 } 



+ 



sm cr„ 



On 



cos a a a b 



L sm a b 



Ob 



Tr 



rpLrpK e -iP a T/2 



2 cosa a Tr 



T L e -iP a T/2 



a K sin a . 



-Tr 



(T 



rpKrpL e -iP a .T/2 



;i70) 



The remaining expression in the traces here are all bounded terms and hence it is sufficient 
to show from the above that none of the terms in the curvature operator diverge. 

One also evaluates the following then using exactly similar techniques, 



Tr 



glgaT^b 



Tr 



9b9'9a 
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